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Abstract 

Evolution process could be calculated from the relativistic hydrodynamic equation with 
certain estimated initial conditions about a single spherical fireball here. So one could estimate 
a kind of initial condition qualitatively with a possible energy density about eo ~ 1.9 GeV / fm 3 , 
based on this process to fit the experimental data at thermal freeze-out. The evolution from 
a cylindrical fireball will be discussed simply in a later chapter. 
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1. Introduction 

PS. 

It has been suggested that in relativistic heavy ion collisions the QGP state would be formed, 
in which quarks and gluons are free to roam within the volume of the fireball created by the 
collisionpp|2][3]0]. Studying the hydrodynamic evolution of the fireball would be helpful to learn 
the QGP and the phase transition from QGP to hadron gas or whether it has happenedjS . The 
evolution of a fireball with certain simple (energy density) distributions in central incident is 
computed to give some coarse estimations of the initial energy density and other information that 
may be the signatures of the QGP[Hj. 

In our calculations, some conditions are presumed. (1) The evolution could be described as a 
quasi-static process, as well as the local thermal equilibrium is formed. Every small mass region 
could be described with the mean energy density, pressure and cm. velocity. (2) Ideal gas is 
provided, therefore p — ^e. (3) The fireball given has a spherical shape and there are no revolving 
fluid movements in central incident, all movements of the relativistic fluid are radial, p, e are also 
radial fields. 

In fact, the original shape is not known exactly especially there is a finite time during the 
collision proceed. The fireball under building also evolves at the same time. For Na49, it is about 
1.5 fm/c. Therefore, spherical model is used here, for the facilities of computing and also due to 
our confidence that the difference could be accepted comparing to uncertainties by other causes in 
the calculations. Cylindrical expansion from a flat fireball will be discussed in Chapter 5. 
With these conditions provided, one can use the relativistic hydrodynamic equation 



where 



(1) 



Tftv = (e + p)u^u v - g^p, 



(2) 



to calculate it. 

2. Hydrodynamic Evolution 

In order to describe the global fireball properly and easily, variables in the form of spherical 
coordinate (r, 9, <p) are used while still keeping the hydrodynamic equation in Minkowski from. 
Because the fluid velocity is always radical as been mentioned, one has (4-velocity) ug = u v = 0, 
and 

u x = u r sin cos (p, 
u y = u r sin 9 sin ip, 
u z = u r cos8. 

Now Eq (1) can be transformed by this way 
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Based on the relations above one gets that, when v = 1, 2, 3, equations are same 



(e + p) 



dv% d(u r Uo) u 



dr 



+ 



dt 



+ 2^ 

r 



,d (e+p) 
dr 



■ u r uo 



d(e+p) dp _ 



dt 



dr 



(4) 



when v = 0, 
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For v = v r = Ur/uo, one has 



and Eqs (7) turn to 
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They are non-linear partial differential equations which could describe the expansion process 
of a spherical fireball. From (9) one can find that the geometrical shape and velocity distribution 
are invariant to the scale of energy density, from the form (<91ne). This is because the equation of 
states of ideal gas (p = -|e) is used. This makes some geometrical data independent of the initial 
energy density, but only depend on its relative distribution (i.e. the value of a) in this model. 



3. Calculations 

Two-step Lax- Wendroff Method |7] is used here to compute the evolution. Gaussian distribution 
and mean distribution of energy density are tried as the initial conditions to run the programmes. 
Gaussian condition works much better than the other in the calculations. 

In order to minimize the computational error and to raise the stability of the programme, a 
small variable e v is added as a correction to the energy density to avoid it too close to zero 
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(10) 



This correction is so small that it has not a little influence on the final results, but it helps the 
programme to compute the evolution for a long time (6-10 fm/c) enough to freeze-out. It could 
be regarded as the energy density of the base state of physical vacuum. 

The programme could output a series of data of energy density, fluid velocity and size. The 
Root Mean Square (RMS) Radius (R rm s) is 
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where E s is the total rest energy 



G 

and Ar is the stepsize about location. The relation to the effective radius (RMS on projection to 
one dimension) is 
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It is interesting that (E s * i^ ras ) looks conservative during the fireball expands, no matter which 
initial conditions it has. It should be a mirror of hydrodynamic equation (energy-momentum 
coservation) . After setting some acceptable conditions on the borders, the programme could evolve 
about 12 fm/c stably. 

4. Experimental Analysis and Evolution Results 

The data can be observed in the experiments are the overall time of expansion t, the transverse 
energy distribution via pseudorapidity dEi/drj, the emission source size a and the transverse 
velocity vt when the the fireball freezes out. Initial inputs are central energy density eo and the 
effective size of its distribution ctq. Trying appropriate eo, <Jq and a fixed t, can fit the other 
experimental data very well. 

The volume of the cross section of central colliding region in cm. frame can be estimated asjS] 

6 sjA\ +A'\ + 2A 1 A' 2 cosh Y L 

where Yj, is the rapidity of incident nucleus in the laboratory frame. 

The radius of this region r c is about 1.5 fm for RHIC and 3.26 fm for SPS 158 GeV 2m Pb. 
According to a Gaussian distribution, this radius should be smaller than the RMS radius and larger 
than the effective radius, that is 

CO < r c < Rrms- 

So from eq (13), for SPS one has 

1.89/m < ct < 3.26/m. 

Initial effective radius will be tried from this domain. 

The overall time of expansion (to freeze-out) t is observed about 8 fm/c near mid rapidity, 
decreasing slightly to 6 fm/c at high rapidity, with a Gaussian radius (mean square error) a = 
Ro = At =3.5 fm/c jHl- We set it at t = 7.5 fm/c as the time when freeze out. 

Transverse energy distribution calculated between 2< 77 <4 shows that there is a peak near 
V _ Vcm. =2.9 and depended on the initial effective size (To, because pseudorapidity is used instead 
of rapidity. The larger the (To is, the larger the central pseudorapidity is and the higher the peak 
is (, about dEx/dr] oc <7o~ 3 m region (1.89/m < <jo < 3.26/m)). At the same time, the initial 
energy density eo contribute the peak value too (, dEx/dr) oc eo). But the location of the peak is 
still invariant to eo, due to our assumption of ideal gas. To fix the location and peak value could 
determine the initial parameters. 

NA49 experiments showedjS] that there is a peak around 77 = 3.0, and the peak value is 405 
GeV per unit, (see Fig 1.) Bjorken formula [TU| gives eo = 3.2 GeV/ fm 3 to fix the peak value. 
Setting eo = 3.2 GeV and cto = 3.13 fm (both are near the upper limits) could also fit the peak 
value well, (see Fig 2,) but the peak location is 77 = 3.133, a bit larger than expected and the peak 
width is only about 1 unit far smaller than the experiment. The value decrease quite rapidly far 
from the centre, while it is still larger than 300 GeV in experiment when 77 = 2 and 4. 

Here cm. velocities of every small regions are used instead of the particle velocities at freeze- 
out to compute the pseudorapidity and transverse energy. Because this method did not mention 
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Figure 1: The Pb+Pb data of NA49. The peak is at 3.0 unit. The peak value is 405 GeV. 



the fluid temperature and its components, one can not calculate the real momentum distribution 
in that small region and give the exact result. For this reason, the real curve should be more 
smooth and the peak should be a little lower. It means the real data are more far from the the 
experiment. Most of all, even if the smooth effect dominate the distribution, the total Et and the 
peak location will not change much. Counted in Fig 2, the total transverse energy between rj = 2.0 
to 4.0, is only about 433 GeV, while in Fig 1, the total transverse energy is no less than 700 GeV. 
Even though set the effective size as 3.2 fm to the limit, the total transverse energy is only about 
464 GeV, while the peak value is 435.4 GeV and its location comes up to 3.147 unit. Larger cto's 
and initial energy densities than those are not suitable. 

So one reasonable explanation is that the experimental datajH] contain a huge background ^1]. 
This background is probably produced by the other small fireballs if multi-fireballs are emerged 
after the collision. The major signal from the central fireball forms the peak and the others make up 
a background. By cutting off the background about 290±20 GeV [H], the peak value decreases to 
about 115±20 GeV and total the transverse energy turns to 123±40 GeV. Choosing a combination 
of ctq and £q properly, could fit it very well. 
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Table 1 Possible eo at a = 2.5 fm 
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Table 2 Possible e at er = 2.6 fm 



Figure 2: The dEx/dr) distribution at eo — 
3.2 GeV/fm 3 and cr = 3.13 fm. The peak 
value is 403.9 GeV and the peak location is 
about 3.13 unit. 



Figure 3: The dEx/dr) distribution at eo = 
1.9 GeV/fm 3 and a Q = 2.55 fm. The peak 
value is 114.0 GeV and the peak location is 
about 3.03 unit. 



Calculations show that from eo near 2.7 GeV/fm 3 at ctq = 2.3 fm (a freeze-out = 4.300 fm), eo 
near 2.0 GeV/fm 3 at <r — 2.5 fm (a freeze-out = 4.350 fm) to eo near 1.8 GeV/fm 3 at cr = 2.6 
fm (ct f r eeze-out = 4.375 fm) are all permissive. Considering that the peak location could not 
be too much larger or smaller than 3.0 and the transverse velocity should not be too high, tr is 
defined to [2.30,2.60]. 

The effective radius at freeze-out in experiment varies from 3.8 fm by proton correlation [T] [T5] 
to 6.5 ± 0.5 fm by pion correlation |8"| |12|. This is quite dramatic and makes it pretty difficult to 
determine the initial size precisely. The effective radius got at freeze-out here varies from 4.15 fm 
to 4.7 fm according to the initial domain ao = 1.9 ~ 3.2 fm (, see Fig 4). This result is a little 
larger than the the data from proton correlations and smaller than the data from pion correlations. 
All could be acceptable, including our conclusion a = 4.3 ~ 4.4 fm which is quite near to the data 
from proton correlations. 



5.50 
5.00 



Z.50 
2.00 



Figure 4: Fireball size during evolution. e =1.9 GeV/fm 3 ; , ctq =2.55 fm. 



Transverse velocity vt= 0.55 is report ed[5| |13j. Again because the local thermal momentum 
distribution and the particle component can not be provided from the hydronynamic method, the 
cm. velocity at the effective radius v a at freeze-out is used to compare with the transverse velocity. 
From (To = 2.3 fm to 2.6 /to, v a reduces from 0.645 to 0.613. We intend to select a relatively 
larger initial effective radius with a lower freeze-out velocity. 



Energy density at freeze-out is estimated about 0.05 GeV/ /m 3 pp. Evolution results with 
parameters discussed above are about 0.051 to 0.065 GeV/fm 3 , very approximate. Smaller and 
larger a will produce too tiny or huge result, although weone can tune the initial energy density 
to give a small correction. 

From above, our selected estimation of the initial parameters is about 1.9 ± 0.3 GeV/ fm 3 at 
(To = 2.55 fm. Detailed results are listed in Fig 5, 6, 7, 8. From Fig 5, one can see that central 
energy density reduces very slow (about 0.1 GeV/ fm 3 ) in the first 1 ~ 2 fm/c. State that the 
possibilities of initial energy from range 1.4 to 2.4 GeV/ fm 3 with relevant effective size could not 
be removed completely either. 
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Figure 5: Energy density evolution at differ- 
ent locations. e = 1.9 GeV/ fm 3 , ctq = 2.55 
fm. 



Figure 6: Energy density distributions at dif- 
ferent time. e — 1.9 GeV/fm 3 , ao = 2.55 
fm. 
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Figure 7: Fluid velocity evolution at dif- 
ferent locations. The velocities are from a 
set locations, but not according to the same 
fluid parts. e = 1.9 GeV/fm 3 , a = 2.55 
fm. 



Figure 8: Fluid velocity distributions at dif- 
ferent time, eo = 1.9 GeV/fm 3 , ao = 2.55 
fm. 



5. Cylindrical Fireball 

In this chapter, a short discuss is done on the evolution of a flat cylindrical fireball. One more 
equation is added to the evolution equations (10) and one dimension of data in memory expand 
too. The hydro dynamic equation(s) could be written to these forms 
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Let t;,. = u r /uo, v z = u z /uq, the partial time forms 
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6. Summary 

Hydrodynamic equation is used to compute a spherical fireball created by the relativistic heavy 
ion collisions. The evolution works very well. It can produce kinds of data to compare with those 
from experiments. While, although these equations do not have any free parameters, but due to 
the complex, unknown and severe uncertain initial conditions, only a qualitative process could be 
given. The estimate of initial data is only a kind of attempt. 

The experimental data is likely to contain a huge background. It is reported that the initial 
energy density could be reduced to 0.91 GeV/ fm 3 \n\, by cutting off the background. To use 
hydrodynamic method to deal with this problem here, the initial energy density is estimated 
about eo ~ 1.9±0.3 GeV/ fm 3 . Thinking that the results are more sensitive to the initial size than 
to the initial energy density, the real error range may be larger. The result is not so striking as 
the estimation of eo = 3.2 GeV/ fm 3 got before. The possibility of the QGP production in CERN 
SPS is still not clear. 
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